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Fig. 2 




CO CO 

Jh > 
CD 

B CD 

CD C/] 



Jh 

CD 



CO 

Jh 
bX) 
O 
Jh 

Oh 












o 






CD 






Jh 


CD 


c 


D 


s 


• I— 1 








CO 


c: 


CD 


Jh 


o 


bjO 


CD 


o 


C 


a 


♦1—1 


CO 


E 






CD 


•T—i 


o 


+-> 


00 







CD 






Jh 


c 






• 1— 1 




at 


CD 




Jh 




CD 


CD 


c 


4-> 


Q. 


CO 


CO 


£ 




CD 


CD 


o jz: 









^ bJO 




CD e 




.2 w 




^ c 


po 


C ° 


Jh 


•rH O 




•i— 1 




CD i — i 


CO 


buO'rH 




£ GO 


co 


CO 


>> 


JT C4H 


Jh 


o o 


O 



measured temperature 





■+-> 




CD 




CO 




CD 




CD bjO 




Jh £ 




=3 CO 




+-> _C 




CO O 


Jh 


Jh 


CD 


CD +-> 




a c 


CO 


E -rH 


0 


CD O 







CD 




Jh 




=3 








ca 




S Jh 




CD CD 




a. 




<n B 




CD 






o 




♦i— i 


r— 1 Jh 


+j 


O CD 


o 


Jh 


CD 


-H> 03 


Jh 


C CD 


Jh 


O _C 


O 


a ' — ' 


O 



Jh 
CD • 



Jh 
O 

CO 

C 

CD 

00 

c 

o 



CD O 

E CD 

CO +-> 

•rH CD 

-a -a 



measured 
diameter 




measured pulling speed 







E 




CD 




+-> 




CO 








CO Jh 


CD 


i ! +-> i 1 


O CD O 


Jh I 


5 Jh 


■P CO H 




o -o o 


o x — ' o 


-o 




CD 


i— i 


CD 


o 


a 


Jh 


co 


+-> 






o 




CJ 


1 — 1 




1 — i 


i — i 




Oh 


a. 




bo 




c 




•1— 1 




1 1 




1 1 




D 








-a 


00 


CD 


a 


CD 


CO 


a. 


CD 


CO 


E 



c 




•iH 


bo 


CD 


c 


tiO 




C 


^H 0 


CO 


^H CD 






o 


a oo 





Jh 


+^ 


CD 


CD 


-Hi 


bXJ 


CD 


Jh 


E 


CO 


CO 


+-> 


•rH 



ML6272 



10/521035 



3/15 
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START 



3 



Calculate temperature 
distribution of silicon single 
crystal rod at time t by 
computer based on 
comprehensive heat 
transmission analysis. 



T 



Input pulling speed 
measured profile of 
silicon single crystal rod 
from start of pulling to 
time t to computer. 



Input set pulling speed of 
silicon single crystal 
rod from start of pulling 
to end of pulling to 
computer. 



Calculate diffusion and 
defect formation of point 
defect by computer using 
defect simulation method 
based on set pulling speed 
and pulling speed 
measured profile. 
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Fig. 5 
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from© 







Decrease set pulling speed by 
predetermined small 
auantitv. 






Calculate diffusion and 
defect formation of point 
defect by computer using 
defect simulation method 
based on corrected 
oulline speed. 
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to® 



Increase set pulling speed by 
predetermined small quantity. 



Calculate diffusion and 
defect formation of point 
defect by computer using 
defect simulation method 
based on corrected 
pulling speed. 






YES 






^YES 


Substitute cor 
speed with set 
and also incre 
temperature. 


rected pulling 
pulling speed 
ase set heater 




Substitute cor] 
speed with set 
and also decre* 
temperature. 


rected pulling 
pulling speed 
ase set heater 
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c 



START 



3 



from® 



Arbitrarily define parameter group (P^ 
?2y • Pn) obtained by varying at fixed 
intervals variable which can be fed back 
to silicon single crystal rod 
manufacturing conditions (pulling speed, 
pulling length Lj and pulling height H]). 



p x =p 



Input coordinate data and physical 
properties of mesh point of each 
member, P x and turbulence parameter 
of silicon melt, to computer. 



Obtain surface temperature 
distribution of each member based 
on calorific value of heater and 
emissivity of each member. 



Obtain internal temperature 
distribution of each member based 
on surface temperature distribution 
and thermal conductivity of each 
member, and then further obtain 
internal temperature distribution 
of silicon melt taking convection 
into consideration. 



Obtain solid-fluid interface shape 
of silicon single crystal rod and 
silicon melt in accordance with 
isothermal line including tri- 
junction of silicon. 
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Obtain coordinate and temperature of 
mesh in silicon single crystal rod and 
input them to competer. 




Divide time from pulling start time t 0 to cooling 
end time t 1 of silicon single crystal rod at 
predetermined intervals At, and obtain pulling 
length L u pulling height ti 1 and internal 
temperature of silicon single crystal rod from 
coordinate and temperature of mesh in silicon 
single crystal rod every At. 



Solve diffusion equation based on diffusion 
coefficients and boundary conditions of vacancies 
and interstitial silicon in silicon single crystal 
rod, and obtain concentration distributions of 
vacancies and interstitial silicon after At. 



Obtain formation start temperature Tvo of void 
based on concentration distribution of vacansies. 

— — j— - 
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to® from® 

| , | 

Obtain formation start temperature T p of high- 
temperature oxygen precipitate based on 
concentration distribuiton of vacancies. 




T : lattice point 

temperature of mesh in 
silicon single crystal 
rod 



Obtain concentration of void when 



Obtain radius of void when T<Tvo. 




YES 



Obtain concentration of high- 
temperature oxygen precipitate when 
T=T P . 



Obtain radius of high-temperature 
oxygen precipitate when T<T P . 
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t 0 = t 0 +A t 




NO 



YES 



Partition P band where high-temperature oxygen 
precipitate exists from non-defective area where 
both void and high-temperature oxygen precipitate 
do not exist, and obtain first isoconcentration 
line HC1 X which partitions non-defective area to 
be positioned on single crystal bottom side apart 
from P band. 







Determine maximum value in first 
isoconcentrtion line HC1 X as Z(HC1 X ). 







Obtain first distribution line BC1 X which in 
isoconcentration line that interstitial silicon 
concentration is 0. 126% with respect to 
interstitial silicon equilibrium concentration 
in silicon melt and partitions non-defective area 
to be positioned on single crystal top side apart 
from B band. 



0. 126% indicates concentration 
for B band observed from 
lifetime map or the like in 
interstitial silicon 
concentration obtained by 
performing simulation. 
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Fig. 10 
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Px=P i+ 1 







Determine minimum value in first 
distribution line BC1 X as Z(BC1 X ). 






Obtain difference A 2 Y (X) between 
Z(HClx) and Z(BC1 X ). 




Display manufacturing condition 
under which AZi(X) becomes 
maximum in group of parameters P^ 
P 2 , ""',Pn in silicon single crystal 
rod manufacturing condition. 



(end) 
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Fig. 11 
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10 silicon single crystal 

pulling apparatus 
13 silicon melt 
24 silicon single crystal rod 
S tri-junction of silicon 
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P band (high- 
temperature oxygen 
precipitate area) 



P band (high- 
temperature oxygen 
precipitate area) 




(V/G*) 

first isoconcentration 
line HC1 X 

first distribution 

line BC1 X 



(V/G*) 

second isoconcentration 
line HC2 X 
second distribution 
line BC2 X 

axial derection of silicon single 
crystal rod( 7 coordinate) 
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Fig. 13 
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